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ABSTRACT 

The system PSR J0737-3039 is the only binary pulsar known to consist of two radio 
pulsars (PSR J0737-3039 A and PSR J0737-3039 B). This unique configuration allows 
measurements of spin orientation for both pulsars: pulsar A's spin is tilted from the 
orbital angular momentum by no more than 14 degrees at 95% confidence; pulsar B's 
by 130 ± 1 degrees at 99.7% confidence. This spin-spin misalignment requires that 
the origin of most of B's present-day spin is connected to the supernova that formed 
pulsar B. Under the simplified assumption of a single, instantaneous kick during the 
supernova, the spin could be thought of as originating from the off-center nature of the 
kick, causing pulsar B to tumble to its misaligned state. With this assumption, and 
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using current constraints on the kick magnitude, we find that pulsar B's instantaneous 
kick must have been displaced from the center of mass of the exploding star by at least 
1 km and probably 5-10 km. Regardless of the details of the kick mechanism and the 
process that produced pulsar B's current spin, the measured spin-spin misalignment 
in the double pulsar system provides an empirical, direct constraint on the angular 
momentum production in this supernova. This constraint can be used to guide core- 
collapse simulations and the quest for understanding the spins and kicks of compact 
objects. 

Subject headings: pulsars: individual (J0737-3039) — supernovae: general 



1. Introduction 



The radio-pulsar syste m PSR J0737-3039 is the only binary pulsar k nown to consist of two radio 
pulsars: PSR J0737-3039 A (jBurgav et alJbooal ) and PSR J0737-3039 B (|Lvne et alJl20o3 >. Tabled] 
gives the parameters of this system. This unique configuration has permitted measurements of spi n 



orientation for both pulsars (jFerdman et al. 



2008 



Lvutikov Thompson 



2005 



Breton et alJl200a) : 



pulsar A's spin i s tilted from the orbi tal angular momentum vector by no more than 14 degrees at 
95% confidence (jFerdman et al J 12008 ): pulsar B's by 130.0l{j de grees at 99.7% confidence. Here 
we argue that this large difference between the two pulsar spin tilts requires that the origin of most 
of B's spin is connected to its supernova (SN) explosion; the spin of B's progenitor, expected to be 
aligned with the pre-SN orbit due to tidal interactions, cannot be invoked to explain the present- 
day misaligned spin of pulsar B. PSR J0737-3039 B is currently believed to have formed from 



an electron-capture supernova triggered in a massive O-Ne-Mg white dwarf (Ivan den H euve] 



Will ems Kalogera 



200£ 



2004; Piran fc Shavh 



2005; 



Stairs et al 



(]; van den Heuvell 120071 : fereton et al. 20081 : Wong et al. 201ol ) . Our results demonstrate that 



2004 



200d : IWang et al.boodJwillems et al 



whatever the details of its formation mechanism, the supernova that formed PSR J0737-3039 B 
produced the majority of its current spin. If the source of the present-day spin of pulsar B is a 
single, impulsive kick, then this kick must be off-center so that it tumbles the pulsar to its current 
orientation. Using constr aints on the SN ki ck magnitude derived from the orbital and kinematic 
parameters of the system (jWong et al.ll2010l ) we find that this kick must have been displaced from 
the center of mass of the exploding star by at least 1 km and probably 5-10 km. Such offset 
distances ar e a significant fraction of t he expected radii of neutron stars. Off-center kicks were first 
suggested in ISpruit Phinneyl (119981 ) on purely theoretical grounds. 



2. Evolutionary History 



PSR J0737-3039 likely evolved from two s tars originally massive enough t o undergo supernova 
(SN) explosions and form two neutron stars (jTauris van den Heuvell 12006) at the end of their 
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nuclear lifetimes. Given the measured spin magnitudes and inferred magnetic fields (IBurgay et al 
20031 : iLvne et al.ll20ol iFerdman etaDbood : iLvutikov & ThompsoiJliooa breton et alibood ). pul- 
sar A was the first-born neutron star, while pulsar B formed in a second SN. After the first SN the 
system passed through a high-mass X-ray binary phase. In this phase, pulsar A accreted matter 
from its companion, leading to some spin-up. Eventually, pulsar A's companion evolved off the 
main sequence and its expanding hydrogen envelope enveloped pulsar A. In this common-envelope 
phase, tidal interactions between the stars circularized the orbit and are expected to have aligned 
the spins of pulsar A and of pulsar B's progenitor with the orbital angular momentum axis (per- 
pendicular to the orbital plane) . The transfer of orbital kinetic energy to the envelope eventually 
removed the outer layers of pulsar B's progenitor, leaving pulsar A in a tight orbit with the ex- 
posed helium-rich core of B's progenitor. After another b rief period of mass transfer onto pulsar A 



(Dcwi h van den Heuvel 



200 



Willems fc Kalogeral 12004 ) . the helium star exploded in the second 



SN, forming pulsar B. As a result of the multiple mass-transfer phases between the two SN events, 
just before pulsar B's SN the system was in a close, circular orbit with both stars' spins aligned 
with the orbital angular momentum vector. 

Due to asymmetries associated with the SN ejecta (matter and/or neutrinos) SNe are thought 
capable o f imparting a signi ficant recoil impulse, a "kick", to any remnant surviving the explosion 



(see, e.g.. iJanka et al.l (|2008l ) and references therein). When a SN occurs in a binary system, these 



kicks can significantly alter the orbital parameters or even disrupt the binary. The kick component 
directed parallel to the pre-SN orbital plane causes a change in the eccentricity and semi- major axis 
of the orbit; the component perpendicular to the pre-SN orbital plane can also cause a change in the 
inclination of the orbital plane. In the PSR J0737-3039 system, puls ar A's small spin-tilt angle (less 
than 14 degrees at 95% co nfidence using a two-pole emission model (jBurgay et al.ll2003l : ILvne et al 



2004; IFerdman et alJl2008f)1 is indic ative of a relatively small out-of-plane kick from the SN that 



formed pulsar B (|Wong et al.ll2010l ). Pulsar A's spin-orbit misalignment occurs only because the 
orbital plane is tilted by the SN kick, while pulsar A's spin remains fixed in the inertial frame 
aligned with the pre-SN orbital angular momentum axis (Figure [1]). Such a spin tilt for pulsar 
A occurs independently of the effects of the second SN on pulsar B's spin. In other words, the 
observed tilt of pulsar A's spin by itself does not require any change in the spin angular momentum 
of pulsar B relative to its progenitor. However, unless the SN contributes significant amounts of 
angular momentum to the nascent pulsar, the orientation of pulsar B's spin will be the same as 
its progenitor's spin, i.e. aligned with the pre-SN orbital plane and pulsar A's spi n. Surprisingly. 



pulsa r B's spin is in fact retrograde: tilted by 130. Ol^ degrees (99.7% confidence: IFerdman et al 



20081 ) relative to the current orbital angular momentum vector (see Figure H]) ! 



3. The Need for Spin Angular Momentum from the Supernova 



To produce pulsar B's retrograde spin, the SN must have significantly torqued pulsar B, caus- 
ing it to tumble to the currently observed spin-orbit orientation. The pre-SN spin, So, the angular 
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momentum produced by the SN ejecta, AS, and the post-SN spir0, Ssn, are related by the con- 
servation of angular momentum 

S 5A r = S + AS. (1) 



To determine AS, we must know Sq and S$Ni but we only know the direction, not the mag- 
nitude, of So and the relationship betwee n Ssn and the spin measured today is complicated by 
relativistic precession ([Breton et al.ll2008l ). However, we can still place constraints on AS. Rela- 
tivistic precession causes the individual pulsar spins to precess about the total angular momentum 
of the system, which is approximately parallel to the orbital angular momentum. Such precession 
preserves the angle between the total angular momentum and the spin (which is the spin colati- 
tude), but not the azimuthal orientation. Thus, the colatitude of Ssn relative to the normal to 
the current orbital plane is equal to the colatitude of the current spin — 130 degrees. Based on 
the spin of pulsar A, the current orbital plane could be tilted at most 14 degrees relative to the 
pre-SN orbital plane. Therefore the colatitude of Ssn relative to the pre-SN orbital plane — and 
therefore relative to So — is at least 116 degrees. This is also the minimum angle between So and 
AS. The angular momentum produced by the SN must be significantly mis-aligned with the pro- 
ge nitor spin. To date, most SN s i mulations have focused on non-rotating progenitor s (for example, 



sec 



Blondin &: Mezzacappal 120071 ; iRantsiou et al.ll201ll ; IWongwathanarat et al.ll2010i ); it remains to 



be seen whether the spin produced by the SN from the collapse of a rotating progenitor can be so 
significantly mis-aligned with the progenitor's rotation axis. 



The typical moment of inertia ( Spruit fc Phinneylll998l ) of a neutron star is 0.36MR 2 ; using 
pulsar B's measured mass of 1.25M Q (see Table [1]) and a radius of 10 km, its current spin angular 
momentum is 2 x 10 45 g cm 2 s _1 . Since this spin is retrograde whereas the pre-SN spin is roughly 
aligned (within 14°, given pulsar A's small spin tilt) with the current orbital plane, we can place a 
lower limit on the change of angular momentum needed to explain pulsar B's large and retrograde 
spin tilt, 

AS > 2 x 10 45 g cm 2 S -\ (2) 

where equality holds when Sq = 0. Because the angle between So and Ssn is greater than 90 
degrees, any progenitor spin only increases the amount of angular momentum that must be added 
to the pulsar by the kick. This is demonstrated geometrically in Figure [2fc. 

The above discussion has been fully general. To extract more constraints from the observed 
spin-spin misalignment, we must make some assumptions about the origin of the pulsar spin. As 
a simplified model to elucidate the scales involved in this scenario, let us assume that the same 
impulsive kick (i.e. linear momentum) that changes the orbit of the system is also offset from the 
center of mass of pulsar B, and therefore applies a torque sufficient to produce the observed spin 



1 It is important to distinguish between pulsar B's spin vector right after the supernova and its pre sent-day spin 



vecto r because relativistic effects cause the spin vector to precess about the total angular momentum ([Breton et al 
2008), leading to a time- varying azimuthal component. 
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angular momentum. The kick and offset vectors must lie in the plane perpendicular to AS (see 
Equation [3j) . The kick velocity, vk, the offset vector relative to the center of mass, r, and the 
change in B's spin vector are related by 

AS = r x Ap = r x M B v K , (3) 

where Ap = Mb^k is the change in linear momentum induced by a change in velocity of vk in an 
object with mass Mb- The offset length r and kick velocity magnitude vk must then satisfy the 
inequality 

M B v K 

where AS is the magnitude of the change of B's spin; the equality holds only when the kick and 
offset are perpendicular to each other. 

The relative orientation of the current spin provides a constraint on the kick direction in this 
scenario. Let the colatitude of AS relative to the pre-SN orbital plane be #a; because the angle 
between So and Ssn is greater than 90 degrees, no matter the magnitude of So we must have 
#A > 116 degrees. Let the plane perpendicular to AS make an angle ip with respect to the pre-SN 
orbital plane. Then we have ip = 180 — #a < 64 degrees. The kick, Vr-, lies in this plane and 
therefore must have colatitude Ok that satisfies 

90 - V = 26 < 9 K < 154 = 90 + V (5) 
This geometry is illustrated in Figure [2j The constraint in Equation [5] is consistent with the 



constraint on kick colatitude in Figure 7 of (jWong et al.ll2010l ). but tighter. 



The constraint on the magnitude of the spin change, Equation ([2]), together with Equation 
@, imply a lower limit on the offset distance 

/25km s^\, 

r > 3.2 km. (6) 

V V K J 

If we assume that the core of pulsar B's progenitor was in synchronous, rigid-body rotation just 
before the supernova then S$ R ~ 2 x 10 45 g cm s , and the limit on the offset distance rises by a 
factor of 1.8: 

/25km s _1 \ 

r S R > 5.8 km. 7 

V V K J 



Wong et al.l (|2010l ) used the measured semi-major axis, eccentricity and proper motion of the 



J073 7-3039 system (see Table [T]) to constrain the kick imparted to the system by the second SN. 



The I Wong et al.l (|2010l ) analysis assumed that the pre-SN orbit was circular and that the system 



came from a progenitor population with number density 



//I R. : ) - /',, cxi") ( -J^j ex P jf 1 ) <*) 
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with h,R = 2.8 kpc and h z = 0.07 kpc the galactic scale length and height, respectively, moving 
with the local galactic rotation velocity. The SN kick and mass loss must then induce the current 
eccentricity and semi-major axis in the orbit and give the system as a whole a velocity such that 
it m oves in the galactic potential to its current location in the 100 to 200 Myr since the second 
SN (jLorimer et al.l 120071 ). Because of the uncertainty in the amount of mass loss, pre-SN semi- 
major axis, pre-SN galactic location, the age of the system, and the measurement uncertainty in 
the current orbi t al pa rameters, a range of kick m agnitudes between and 60 km/s is allowed in 
the lWong et al.l (|2010i ) analysis (jWong et al J 120101 . Figure 5). In Figure [3] we show the probability 
distribution of minimum offset distances implied by this distribution of kick velocities. Even for 
large kick velocities, the minimum offset distance exceeds ~ 1 km. For the smallest allowed kicks, 
the minimum offset distances exceed the current ~ 10 km radius of the neutron star. 

The magnitude of the offset, r, required to produce the needed AS depends on the relative 
orientation of the offset and kick, v%. If the angle between r and in the plane perpendicular 
to AS is 6 r K, then 

A .<? 

(9) 



Mbvk sin 9 r K ' 

which is larger than the minimum offset distance by a factor of (sin^ r ^-) _1 (see Equation [3]). 
Kicks that are nearly aligned with the radial vector (small 6 t k) require arbitrarily large off-center 
distances to match the current spin orientation. Even modest misalignments of 40 degrees give an 
enhancement factor of sin -1 40 ~ 1.6 over the offset for perpendicular r and vk- 



4. Discussion 



Multi-dimensional simulations (jBlondin fc Mezzacappal 120071 ) of core-collapse SNe have shown 



that an instability in a stationary accretion shock (SASI) may provide a method of depositing a 
substantial amount of spin angular momentum (2 x 10 47 g cm 2 s _1 ) onto a proto-neutron star as 
part of the co llapse process itself and separa te from any rotation of the progenitor. (In fact, the 



simulations of lBlondin fc Mezzacappal ([20071 ) used non-rotating progenitors.) This would be more 



B (2 


x 10 45 g cm 2 s 1 ). The general consensus from recent evolutionary studies ( 


van den Heuvel 


2004 




Willems & Kalogera 


2004 




Podsiac 


lowski et al. 


2005; 


Stairs et al. < 


>006; 


Wang et al. 


2006; 


Willems et al. 


2006 




van den Heuvel 


2007 




Breton et al. 


2008 


; 


Wong et al. 


2010 


) is that pulsar B 



was formed in an electron-capture SN, where an iron core is never formed and instead core collapse 
is initiated through electr o n captures onto Ne/Mg nuc l ei (IDessart et al.ll2006l ; iKitaura et al.ll2006l ; 



Podsiadlowski et all 12004 : iMiyaji et al. 



1980 



Nomotd 1 1984 1 19871 ) . During the course of such a 
100 km to ~ 10 km. 



collapse, the proto-neutron star shrinks from a radius of - 

We emphasize that the assumpti on of the foregoing discu ssion that the kick is applied at a 
single location is simplistic (see, e.g., ISpruit k, Phinnevl (|1998l )); in a real supernova, both linear 
and angular momentum will be accumulated by the proto-neutron star throughout its formation 
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(a) Pre-SN orbit (b) Post-SN orbit with an on-centre kick (c) Observed Post-SN orbit 




Fig. 1. — Effect of SN kick on binary orbit. The pre-SN orbit containing pulsar A and pulsar B's 
progenitor is shown in (a). The effect of an on-center SN kick that slightly changes the inclination 
of the orbit is illustrated in (b). Notice the post-SN alignment of the two pulsars' spin axes. Part 
(c) illustrates the present-day orbit with a 130 degree misalignment between pulsar B's spin axis 
and the orbital axis. 



Table 1 : J07373-3039 Sys tem Parameters. Except as noted, properties are given in I Bur gay et al 
(|2003h:lLvne et al.l (fcooj). 



Distance 


600 pc 


Galactic Latitude 


245.2 deg 


Proper Motion 


10 km/s 


Spin Period (A) 


22.7 ms 


Spin Period (B) 


2.8 s 


Mass (A) 


1.34 M 


Mass (B) 


1.25 M 


Spin-orbit misalignment (A) 


< 14 deg (95% confidence) (Ferdman et al. 2008) 


Spin-orbit misalignment (B) 


130.0ljj deg (99.7% confidence) 

(Xvutikov & ThomDSon 2005: Breton et al. 2008) 


Orbital Period 


2.4 hrs 


Semi-major Axis 


1.26 R® 


Eccentricity 


0.0878 
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Fig. 2. — Geometry of the angular momentum change due to an off-center kick imparted to the 
nascent pulsar. The x-y-z reference frame is anchored on pulsar B at the time of its SN; the x-y 
plane is the pre-SN orbital plane. In panel (a) we show the relationship between So, Ssn, and AS 
and the orientation of the plane orthogonal to AS (which is inclined by an angle ip with respect 
to the orbital plane). In panel (b) we show that vectors lying in the plane orthogonal to AS — 
like vk — can have an inclination with respect to the orbital plane that varies between — ip and ip, 
leading to the constraint on the colatitude of vk of 90 — ip < k < 90 + ip found in Equation ([5]) . In 
panel (c) we show that as the magnitude of So increases, AS increases in magnitude (see Equation 
([2])) and tilts toward the south pole, reducing tp. Therefore, the most conservative constraints on 
Ok are obtained when Sq = 0, giving 26 < Ok < 154 as in Equation ([5]). 
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10 



r (km) 



Fig. 3. — Distribu tions of minimum offset distances implied by the current orbital constraints 



(|Wong et al.ll2010l ) on kick magnitudes (see Equation ([6])). The solid line gives the distribution 
assuming that So = 0, with a kick that is orthogonal to the offset vector. The dashed line assumes 
that the core of pulsar B's progenitor was in synchronous, rigid-body rotation just prior to the 
supernova (Sq = Sq R ~ 2 x 10 45 gem 2 s _1 ) and that the kick-offset angle is 40 degrees. 
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at many different locations. In this more realistic context, the constraints above on offset distances 
should be interpreted instead as constraints on the offset scale at which the bulk of the linear 
and angular momentum is accumulated. More precise constraints will require detailed modeling of 
the hydrodynamic process of momentum accumulation in the supernova that formed PSR J0737- 
3039B. Nevertheless, it is interesting that the location of kicks inferred from such a simple model 
is consistent with kick origins in the bulk of the shrinking proto-neutron star during the supernova 
(see Figure [3]). Some recent SN modeling suggests that the processes that produce the kick and 
those that impart rotation to the resulting neutron star produce independent kicks and spins, 
and therefore there is little correlation bet ween the kick magnitude and direction and the rota tion 

In 



imparted to the post-SN compact object (jWongwathanarat et al. 



2010 



Rantsiou et al 



2011 



this case the offset length scale inferred above from the dynamical constraints on the kick would 
not be relevant. We conclude that only if pulsar B's spin is actually linked to the torque induced 
by the physical mechanism producing the kick it must be offset from the center of mass of the 
collapsing neutron star progenitor. 

Regardless of the specifics of the collapse process, however, the expected alignment of the spin 
of pulsar B's SN progenitor with the pre-SN orbital angular momentum and the observed mis- 
alignment of pulsar B's spin and orbit at present uniquely imply that pulsar B's spin is dominated 
by angular momentum produced during the SN process, not angular momentum provided by the 
progenitor. The realization of this empirical constraint on angular momentum production in su- 
pernova e presented here i s uniquely enabled by the spin spin misalignment in the PSR J0737-3039 
system (ILyne et al.l I2004J : iFerdman et al.1 120081 : iLyutikov Thompson! 120051 ) and can be used to 
guide core-collapse simulations and the quest for the understanding of compact object spins and 
kicks. 
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